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Abstract — Waveform-diverse  radar  arrays  have  been  proposed 
as  a  method  to  facilitate  single  pulse  imaging  as  well  as  to 
potentially  enable  simultaneous  multi-mode  operation. 
Transmitting  different  waveforms  on  the  elements  of  a  uniform 
linear  array  consequently  raises  the  spatial  and  temporal 
sidelobes  of  the  receiver  matched  filter.  In  this  paper  a 
recursive  minimum  mean  square  error  based  receiver  design 
denoted  as  Space-Range  Adaptive  Processing  is  presented.  The 
new  method  is  capable  of  mitigating  space  range  sidelobes 
thereby  providing  enhanced  sensitivity  for  this  transmission 
scheme. 


L  Introduction 

The  waveform-diverse  array  is  a  conceptual  method  for 
spatially  distributing  energy  from  a  radar  antenna  array  [1-5]. 
The  spatial  distribution  of  energy  raises  the  spatial  sidelobes 
of  the  matched  filter  receiver  due  to  the  increased  energy 
transmitted  in  the  off-“mainbeam”  directions  (relative  to  the 
standard  approach  of  transmitting  an  identical  waveform  on 
each  element).  However,  the  distributed  energy  allows  the 
radar  to  simultaneously  image/search  the  illuminated  scene  as 
opposed  to  electronically  scanning  over  the  same  area.  In 
addition,  the  added  diversity  allows  the  radar  to  perform 
different  functions  in  desired  directions  e.g,  emitting 
communications  and  electronic  attack  waveforms  while 
transmitting  a  sensing  waveform  in  other  directions.  Recently, 
a  number  of  approaches  have  been  suggested  for  waveform 
diverse  array  receiver  designs  within  the  framework  of  MIMO 
radar  [6-8]. 

In  this  paper  a  new  minimum  mean  squared  error  (MMSE) 
based  approach  entitled  Space-Range  Adaptive  Processing 
(SRAP)  is  applied  to  the  waveform  diverse  array  scenario  in 
an  effort  to  jointly  suppress  the  range  and  spatial  receive 
sidelobes  thus  providing  enhanced  sensitivity  for  this 
transmission  scheme.  SRAP  can  be  viewed  as  a  multi¬ 
dimensional  extension  of  Adaptive  Pulse  Compression  (APC) 
[9]  and  Re-Iterative  Super  Resolution  (RISR)  [10]  that 
addresses  the  issue  of  adapting  in  both  the  range  and  spatial 
dimensions.  In  addition,  multi-dimensional  adaptivity 
facilitates  the  realization  of  space-range  nulls  analogous  to 
space-Doppler  nulls  demonstrated  with  Space  Time  Adaptive 
Processing  (STAP)  [11].  Although  several  similarities  can  be 


made  to  STAP,  the  SRAP  algorithm  utilizes  a  structured 
covariance  matrix  to  perform  scene  estimation  where  as  STAP 
typically  relies  on  a  sample  covariance  matrix  formed  from 
measured  data  to  achieve  adequate  clutter  cancellation  for 
moving  target  indication. 

SRAP  is  an  iterative  approach  that  requires  an  a  priori 
estimate  of  the  illuminated  scene  {e.g.  via  matched  filtering). 
The  previous  estimate  is  used  to  construct  a  specific  filter  for 
each  range-angle  cell  of  interest.  The  set  of  filters  are  then 
applied  to  yield  an  improved  estimate  of  the  illuminated  scene 
which  can  subsequently  be  used  to  create  a  new  set  of  filters. 
The  algorithm  has  been  observed  to  converge  after  three  or 
four  adaptive  stages. 


II.  Space-Range  Adaptive  Processing 


A.  Signal  model 

The  waveforms  transmitted  from  an  M  element  uniform 
linear  array  can  be  represented  as  the  NxM  matrix  S  of 

which  the  column  contains  the  length-A  discretized 

waveform  transmitted  from  the  element  of  the  array. 
Note  that  for  traditional  beamforming  the  columns  of  S 
would  be  identical  aside  from  the  elemental  phase  shift  used 
for  transmit  beam-steering. 

The  discretized  model  for  the  received  signal  from  the 
range  cell  on  the  M  elements  of  a  waveform-diverse  uniform 
linear  array  can  be  denoted  as  the  IxM  vector 
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of  ^  element  spacing,  and  z(f)  is  a  IxM  vector  of 
additive  noise  samples. 

The  collection  of  N  temporal  snapshots  of  (1)  can  be 
expressed  as 

Y(f)=  +Z(^),  (2) 

_  9 
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contains  the  complex  amplitudes  of  scatterers  within  2N -\ 
range  cells  of  x{^£,9) .  The  final  space-range  signal  model  is  a 

reorganized  version  of  (2)  and  is  expressed  as  the  NMx\ 
vector 

_  d 

where  (8)  denotes  the  Kronecker  product  and 
z(f)  =  vec^Z^  (f)j .  Using  the  formulation  in  (4),  the 
application  of  a  space-range  matched  filter  takes  the  form 

=  y(^)>  (5) 

in  which  (•)^  denotes  the  complex-conjugate  transpose  (or 
Hermitian)  operation.  Likewise,  the  SRAP  estimate  can  be 
obtained  by  applying  the  adaptive  filter  w  (f ,  ^)  as 

•«SRAp(A6')  =  w(f,6»)^y(f),  (6) 

which  is  derived  in  the  following  section. 

B.  Space-Range  Adaptive  Processing 

The  MMSE  cost  function  for  the  complex  amplitude  in  the 
range-angle  cell  corresponding  to  delay  i  and  spatial  angle  6 
is  given  as 

J{£,e)  =  E  \x{£,e)-yv^  {£,e)y{£)^  ,  (7) 

where  E\^\  is  the  expectation  operator  and  w(f,^)  is  the 
adaptive  filter  for  the  {^,0)  range-angle  cell.  Minimization  of 
(7)  with  respect  to  w*(f,^)  yields 

ys{£,e)  =  (£[y(f)y''  {e)j' e[x*  {£,9)y{£)],  (8) 

in  which  (•)  denotes  complex  conjugation. 


Assuming  the  range-angle  cells  are  uncorrelated  with  each 
other  and  with  the  noise,  the  filter  in  (8)  can  be  expressed  as 

f 

w{£,9)  =  p{£,0)  y)(R(f,(Z>)<H)v^v^)-l-Rz(f)  Sv^Ov^, 

W  ) 

(9) 

where  p{l,0)  =  E  |x(f,^)|^  is  the  power  in  the  range-angle 

cell  corresponding  to  delay  i  and  spatial  angle  6 , 
Rz  (f )  =  (yz^NMxNM  is  noise  covariance  matrix  under  the 
assumption  of  white  noise  (where  (7^  is  the  noise  power),  and 

R{£,<P)=  X  +  (10) 

n=-N+l 
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n>t) 

consists  of  shifted  versions  of  the  composite  waveform 
transmitted  in  the  (j)  direction  according  to  Sv^  .  Note  that 

due  to  the  inherent  transmit  coupling  of  space  and  range,  the 
composite  waveform  is  different  for  different  spatial  transmit 
directions.  Additionally  a  normalization  procedure,  similar  to 
that  used  in  [12],  is  applied  to  the  filter  in  (9)  yielding 


i£,e)  =  - 


w(f,6>)  I  ( 

'(£,9)(Sy,®y,)J(Sy,Oy,f 


in  which  G  =  max[(Sv^(8)v^)^(Sv^(8)v^)].  This 

normalization  alleviates  gain  variations  among  the  set  of 
range-angle  filters  thus  resulting  in  improved  performance 
and  robustness  to  high  dynamic  range  scenarios. 


III.  Implementation 

SRAP  utilizes  a  recursive  MMSE  structure  to  alternate 
between  estimating  range-angle  specific  filters  and  the  range- 
angle  scattering  coefficients  of  the  illuminated  area.  A  block 
diagram  depicting  the  implementation  of  the  SRAP  algorithm 
is  shown  in  Fig.  1.  The  matched  filter  from  (5)  can  be  used  to 
obtain  an  initial  estimate  of  the  scattering  coefficients  which 
are  then  used  to  construct  the  covariance  matrix  in  (10)  for 
each  cell  needed  to  form  the  corresponding  adaptive  filter. 
The  new  estimates  obtained  by  applying  these  unique  range- 
angle  filters  can  be  utilized  to  update  the  signal  covariance 
matrix  after  which  a  new  set  of  filters  can  be  computed  and 
applied.  The  algorithm  generally  converges  after  three  or  four 
adaptive  stages. 
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Figure  1.  Block  diagram  of  SRAP  implementation 

IV.  Transmission  scheme 

Several  approaches  have  been  examined  to  optimally 
achieve  a  desired  transmit  beampattem  via  waveform 
diversity  [13-15].  In  this  paper  we  employ  a  simple  yet 
effective  beam  spoiling  technique.  However  it  should  be 
noted  that  improvements  in  the  characteristics  of  the 
transmitted  waveforms,  i.e.,  low  range/spatial  sidelobes,  will 
yield  enhanced  sensitivity  in  the  receiver  output  for  both  the 
matched  filter  and  SRAP.  Here,  waveform  multiplexing  and 
an  intra-pulse  steering  will  be  utilized  to  create  a  spoiled 
transmit  beampattem.  This  technique  allows  the  user  to 

define  a  set  of  K  angle-specific  waveforms  where  the 
length-  waveform  is  denoted  and  6j^  is  the 

associated  angle.  The  initial  step  in  this  particular  scheme  is 
to  interleave  the  desired  waveform  set  in  a  time-division 
multiple-access  (TDMA)  format  [16]  and  repeat  the  sequence 
to  form  a  matrix  in  which  each  column  is  associated  with  the 
transmitted  waveform  on  the  corresponding  antenna  element 
denoted  as 


r 

Ar 

V 

jEi. 

K  y 

Next,  intra-pulse  steering  is  applied  to  which  imposes 

the  corresponding  elemental  phase-shift  associated  with  each 
of  the  multiplexed  waveforms  denoted  as 
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Traditionally  a  linear  phase  shift  is  induced  between  array 
elements  over  an  entire  pulse  to  transmit  energy  in  a  single 
direction  whereas  in  (14)  the  elemental  phase  shift  is  altered 
multiple  times  during  the  pulse  to  steer  different  portion  of  the 
temporal  waveform  in  different  spatial  directions.  Note  that 
each  discrete  sample  of  an  underlying  waveform  is 
transmitted  in  the  same  direction.  The  spatio-temporal 
structure  of  the  transmitted  signals  may  potentially  enhance 
cooperative  passive  bistatic  radar  while  complicating  passive 
exploitation. 

The  transmit  strategy  can  be  characterized  in  terms  of  an 
aggregate  beampattem  [13]  denoted  as 

5(^)  =  vXSv,,  (15) 

which  can  be  decomposed  into  the  contributions  from  each 
temporal  sample  of  S  resulting  in  the  instantaneous  time- 
varying  beampattem 

B^y{n,d)  =  y'lJXyg,  (16) 

where  is  the  column  of  .  Note  that  the  stmcture 
employed  in  (13)  and  (14)  can  be  augmented  to  create  several 
different  beampatterns. 

As  an  example,  consider  a  M  =  20  element  uniform 
linear  array  with  K  =  \0  waveforms  of  length- =  20 

(thus  10  different  spatial  beams  will  be  formed  during  the 
pulse  interval).  In  this  example  the  waveforms  are  unique 

random  polyphase  codes  and  = -30” for 

k  =  0,-- ,K-\ .  The  resulting  aggregate  and  time-varying 
beampatterns  are  shown  in  Figs.  2  and  3,  respectively.  Both 
beampatterns  are  peak-normalized  to  unity. 


Figure  2.  Aggregate  beampattem  for  waveform  diverse  array 
example 
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.  Time-varying  transmit  beampattern  in  dB  for 
waveform  diverse  array  example 


V.  Simulation  Results 

Again  consider  an  M  =  20  element  uniform  linear  array 
with  half-wavelength  spacing.  For  this  case  the  transmitted 
waveforms  are  those  described  in  the  previous  section.  First 
an  imaging  scenario  will  be  inspected  and  then  a  surveillance 
setting  containing  multiple  targets  distributed  in  range  and 
angle  will  be  examined.  In  all  cases  the  scatterers  are  modeled 
as  point  targets  with  random  phase  and  amplitude  in  additive 
white  Gaussian  noise.  The  scenes  are  constructed  using  181 
spatial  bins  distributed  evenly  between  ±90°  relative  to 
boresight.  The  matched  filter  and  SRAP  will  utilize  9 1  spatial 
bins  evenly  distributed  between  ±90°  relative  to  boresight  for 
processing.  Note  that  the  illuminated  scenes  will  contain 
targets  that  lie  between  the  spatial  directions  considered  for 
processing  resulting  in  spatial  steering  vector  mismatch.  A 
post-processing  equalization  is  applied  to  both  the  matched 


filter  and  SRAP  estimates  to  remove  the  weighting  imposed 
by  the  transmit  beampattern.  Note  that  all  images  are 
normalized  by  the  processing  gain  and  presented  in  dB. 

Figure  4  displays  the  ground  truth  for  the  imaging  case 
which  consists  of  several  closely  spaced  shapes.  Before 
processing  the  SNR  of  the  shapes  range  from  -12  dB  to  0  dB 
(the  total  processing  gain  is  26  dB).  The  matched  filter  output 
in  Fig.  5  is  plagued  with  space-range  sidelobes  thereby 
resulting  in  a  very  poor  image.  In  contrast,  Fig.  6  displays  the 
SRAP  output  which  has  suppressed  nearly  all  of  the  space- 
range  sidelobes  to  the  level  of  the  noise.  A  total  of  four 
adaptive  stages  were  utilized  to  produce  the  image  in  Fig.  6. 

The  "one  size  fits  all"  approach  of  deterministic  filtering 
attempts  to  construct  a  filter  that  is  appropriate  for  all 
scenarios,  typically  by  minimizing  the  peak  or  average 
sidelobe  level.  In  contrast,  the  adaptive  filter  structure  allows 
a  unique  filter  to  be  designed  for  each  range-angle  cell  based 
on  an  estimate  of  the  surrounding  scatterers.  This  approach 
alleviates  the  constraints  on  the  filter  by  requiring  only 
enhanced  performance  for  the  particular  scenario 
corresponding  to  each  filter's  location.  For  example.  Fig.  7 
displays  the  peak-normalized  range-angle  response  for  the 
SRAP  filter  from  range  cell  215  and  angle  0°  after  the  final 
stage.  The  filter  response  indicates  the  output  of  a  specific 
filter  when  multiplied  by  delayed  versions  of  the  spatial 
steering  vectors.  Notice  the  space-range  nulls  at  the  range 
and  angle  offsets  corresponding  to  the  surrounding  scatterers 
thus  demonstrating  the  ability  of  SRAP  to  manipulate  the 
available  degrees  of  freedom  to  reduce  contributions  from 
range/spatial  sidelobes  where  appropriate.  Note  that  the  filter 
response  still  contains  sidelobe  levels  near  that  of  the 
matched  filter  response.  However,  the  SRAP  filter's 
sidelobes  reside  in  locations  that  correspond  to  regions  of  the 
scene  that  do  not  contain  any  large  scatterers  and  thus  do  not 
impact  performance. 
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Figure  4.  Ground  truth  for  imaging  scenario 
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Figure  5.  Matched  filter  estimate  for  imaging  scenario 
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Figure  6.  SRAP  estimate  for  imaging  scenario 
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Figure  7.  SRAP  filter  response  for  imaging  scenario 
(constructed  with  the  adaptive  filter  from  ^  =  215,  ^  =  0°) 


The  radar  parameters  for  the  second  scenario  are  the  same 
as  before  but  the  illuminated  scene  now  contains  9  distributed 
targets  between  range  cells  125  and  175  and  angles  -30°  and 
+30°  as  described  in  Table  1  (stated  SNR  values  are  before 
processing).  After  processing  the  SNR  of  the  largest  target  is 
35  dB.  Similar  to  the  previous  case  the  matched  filter  results 
in  limited  sensitivity  as  illustrated  in  Fig.  8,  in  which  it  is 
difficult  to  separate  the  targets,  whose  locations  are  denoted 
by  the  white  circles,  from  the  sidelobes.  Figure  9  displays  the 
result  after  three  adaptive  stages  of  SRAP  resulting  in 
complete  sidelobe  mitigation  thus  revealing  all  9  targets.  The 
output  power  levels  in  Figs.  8  and  9  are  normalized  by  the 
processing  gain  such  that  the  output  noise  power  is 
approximately  0  dB.  This  particular  mode  enables  wide  area 
coverage  for  a  sustained  period  of  time  thus  allowing  for 
inverse  synthetic  aperture  imaging  to  be  accomplished  at  the 
same  time  as  detection  and  tracking  operations.  It  should  be 


noted  that  different  types  of  waveforms  (e.g.  track,  search, 
imaging)  can  be  transmitted  simultaneously  using  the 
transmission  scheme  presented  here. 

Table  1 .  Target  Description  for  Surveillance 
Scenario 


Range  Cell 

Angle  (Deg) 

SNR  (dB) 

126 

-21 

-10 

129 

0 

9 

137 

-26 

-9 

I4I 

-3 

-11 

144 

20 

7 

146 

26 

-11 

148 

13 

-8 

157 

-10 

-10 

I6I 

-20 

4 
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The  closely  related  one-dimensional  MMSE  spatial 
estimation  technique  RISK  is  capable  of  achieving  finer 
resolution  than  the  matched  filter,  thus  we  expect  to  see 
similar  performance  from  the  coupled  domain  processor 
SRAP.  After  examining  the  spatial  extent  of  the  targets  in 
Figs.  7  and  8  it  appears  that  SRAP  achieves  a  finer  angular 
resolution  than  the  matched  filter,  however  full  evaluation  on 
the  super-resolution  capabilities  of  the  algorithm  is  left  as 
future  work.  Additionally,  in  Fig.  7  some  of  the  targets  do  not 
fall  on  the  spatial  sampling  grid  used  for  processing  (targets 
with  odd  angle  values  in  Table  1)  resulting  in  spatial  steering 
vector  mismatch.  However,  the  SRAP  algorithm  is  still 
effective  at  suppressing  the  sidelobes  of  these  targets. 


VI.  Conclusions 

A  new  minimum  mean  squared  error  (MMSF)  based 
technique  denoted  as  Space-Range  Adaptive  Processing 
(SRAP)  has  been  proposed  that  is  capable  of  mitigating  the 
joint  space-range  sidelobes  inherent  to  waveform-diverse 
arrays.  In  the  spirit  of  the  well  known  STAP  formulation, 
SRAP  utilizes  a  range-angle  coupled  signal  model  allowing 
for  simultaneous  adaptation  in  the  spatial  and  range 
dimensions.  SRAP  has  been  shown  to  exhibit  enhanced 
sensitivity  when  compared  to  the  matched  filter.  In 
combination  with  the  transmission  scheme  presented  here, 
SRAP  facilitates  the  potential  realization  of  some  forms  of 
simultaneous  multi-mode  operation. 
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